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From enemies to friends: chiral symmetry on the lattice 
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The physics of strong interactions is invariant under the exchange of left-handed and right- 
handed quarks, at least in the massless limit. This invariance is reflected in the chiral symmetry 
of quantum chromodynamics. Surprisingly, it has become clear only recently how to implement 
this important symmetry in lattice formulations of quantum field theories. We will discuss real- 
izations of exact lattice chiral symmetry and give an example of the computation of a physical 
observable in quantum chromodynamics where chiral symmetry is important. This calculation 
is performed by relying on finite size scaling methods as predicted by chiral perturbation theory. 

1 Introduction 

Nature as well as physicists like symmetries. The -somewhat amusing- reason for this is 
that symmetries can be broken. A very important concept is the spontaneous breakdown of 
a symmetry: here some symmetry is broken down to a situation with less symmetry when 
tuning a parameter (e.g. the temperature or the coupling strength) of the theory to a critical 
value. In the process a non-vanishing vacuum expectation value is developed that breaks 
the symmetry of the interaction. The process is accompanied by the appearance of so- 
called Goldstone particles (spinwaves) that are massless. This phenomenon comes under 
the name of the Goldstone theorem. An example is the spontaneous magnetization: a metal 
at high temperature is in a symmetric state - the elementary magnets or spins can point in 
any direction such that the net magnetization vanishes. Decreasing the temperature below 
some critical value, a spontaneous magnetization occurs, the spins point into a prefered 
direction and the metal becomes magnetic. 

Within the standard model of elementary particle interactions we know two places 
where such a spontaneous symmetry breaking is supposed to have occurred: in the elec- 
troweak sector spontaneous symmetry breaking manifests itself in the Higgs phenomenon 
with the development of a Higgs field expectation value ($), giving mass to the elementary 
particles, and the appearance of Goldstone particles leading to the W- and Z-bosons. 



In our theory of strong interactions, quantum chromodynamics (QCD), it is a chiral 
symmetry that is assumed to be spontaneously broken. This symmetry allows for an in- 
terchange of left handed and right handed quarks while leaving physics invariant - at least 
when these quarks are massless. In this case a scalar quark-antiquark qq condensate is 
developed and the Goldstone particles are identified with the light pions that are observed 
in nature. 

As stated above the occurrence of spontaneous symmetry breaking is an assumption. 
The phenomenon is inherently non-perturbative and cannot be addressed with approxima- 
tive methods like perturbation theory. However, even with numerical simulations it is dif- 
ficult to test, whether a certain model exhibits spontaneous symmetry breaking (SSB). The 
reason for this becomes clear when the way to detect SSB is considered. Let us choose a 
system that has a finite physical volume V as would be required for numerical simulations. 
Further, we couple the system to an external magnetic field. Spontaneous symmetry break- 
ing is tested in a double limit, where first the volume of the system is sent to infinity and 
then the external magnetic field is sent to zero. If a non-vanishing magnetization remains, 
spontaneous symmetry breaking is identified. Obviously, such a procedure is unfeasible 
within the approach of numerical simulations. 

The way out is the use of chiral perturbation theory 1 . In this approach chiral symme- 
try breaking is taken as an assumption with the consequences of the appearance of non- 
vanishing field expectation values and Goldstone particles. A special situation arises when 
the size of the box becomes comparable to or even smaller than the Compton-wavelength 
of the Goldstone particle. Then the corresponding field can be considered as being uni- 
form and it is possible to set up a systematic expansion that starts in the lowest order with 
an effective lagrangian of this constant mode and then taking systematically higher order 
fluctuations into account 2 . 

2 Example of $ 4 -theory 

Let us give an example of the Ginsburg-Landau or $ 4 theory with 0(A)-symmetry in four 
dimensions. The action of this theory is defined by 



with $ a A-component vector, m the bare mass, A the bare quartic coupling and j a a 
constant external source in direction a of the group O(N). 

Expectation values of observables are computed through the partition function or path 
integral in finite volume V = L 4 , 



Symmetry breaking is detected via a spontaneous magnetization (<&) ^ in the double 
limit 
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Chiral perturbation theory can provide now a prediction for the behaviour of ($)j,y, i.e. 
the expectation value in finite volume and at non vanishing external field. In a situation 
when the external field j and 1/ V are very small, (<&) j,y is given in terms of the expectation 
value of (<j>) at j = and V = oo: 



v = f(u) (4) 

where u is a scaling variable 

u = ($)jV (5) 

that contains the infinite volume and j = spontaneous magnetization (<j>). 
Let us give for completeness the function f(u): 



/(u) = with jy(ti) = - — 

J V U il 



, . (6) 

i(«) 

with Ii, h modified Bessel functions. The important point here is that these functions only 
depend on a single scaling variable u and hence on the quantity of interest, the magnetiza- 
tion ($). 

In an already quite old work 3 the prediction of chiral perturbation theory, eq. (4), was 
confronted with numerical data. In this test the data obtained through numerical simula- 
tions were described very well by the theoretical formulae of chiral perturbation theory 
and the method of chiral perturbation theory turned out to be very fruitful, at least in this 
example of a relatively simple model. The lesson that is to be learnt from this example 
is that finite size effects can actually be used to determine properties of an infinitely large 
system. In this sense, the finite volume system can be regarded as a probe of the target the- 
ory in infinite volume. This gives an entirely new perspective on finite size effects: instead 
of being afraid of them, they can be used to determine important physical information. 



3 The example of quantum chromodynamics 

Chiral perturbation theory as well as the lattice method was developed for quantum chro- 
modynamics in order to deepen our understanding of the strong interactions. Still, until 
relatively recently, both approaches could not really come together, at least in the region 
of very small quark masses where chiral symmetry starts to get restored. The reason was 
that the lattice seemed to be lacking the concept of chiral symmetry and for many years the 
infamous Nielson-Ninomiya theorem 4 was telling us that it would even be impossible to 
implement chiral symmetry in a consistent way on a lattice. 

The situation only changed a few years ago, when an old work by Ginsparg and Wil- 
son 5 was rediscovered 6 . The Ginsparg- Wilson paper contained actually a clue for an- 
swering the problem of chiral fermions on the lattice. The interaction of the fermions is 
described by some particular operator, the Dirac operator, the details of which should not 
be discussed here. Now, in the continuum theory this Dirac operator anticommutes with 
a certain 4-dimensional matrix 75 = diag(l, 1, — 1, —1). On a lattice with non-vanishing 
lattice spacing a, such an anti-commutation property cannot be demanded. If the anti- 
commutation property is insisted on, the fermion spectrum of the lattice theory does not 
correspond to the one of the target continuum theory. 
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The suggestion of Ginsparg and Wilson was to replace the anticommutation condi- 
tion by a relation (now known as the Ginsparg-Wilson (GW) relation) for a lattice Dirac 
operator D: 

l5 D + D lb = aD l5 D . (7) 

Clearly, in the limit that the lattice spacing vanishes the usual anti-commutation relation of 
the continuum theory is recovered. 

The fact that renders the relation eq. (7) conceptually extremely fruitful is that it implies 
an exact chiral symmetry on the lattice even if the value of the lattice spacing does not 
vanish 7 . The notion of a chiral symmetry on the lattice is a conceptual breakthrough and 
renders the lattice theory in many respects to behave like its continuum counterpart with 
far reaching consequences. 

However, as nice as the theoretical progress that followed the rediscovery of the 
Ginsparg-Wilson relation was as much of a challenge are realizations of operators D that 
satisfy the Ginsparg-Wilson relation (see 8 ' 9 ' 10 for reviews). Let us give a particular exam- 
ple for such a solution as found by H. Neuberger 11 from the overlap formalism 12 , based 
on the pioneering work of D. Kaplan 13 . To this end, we first consider the standard Wilson 
Dirac operator on the lattice: 



Av^K^ + V^-aV^} (8) 

with V p , V* the lattice forward, backward derivatives, i.e. nearest neighbour differences, 
acting on a field 

V M $(x) = + n) - $(x) 
V*$(x) = $(x) - $(x + fi) . 

We then define 

A = 1 + s - D w (9) 

with < s < 1 a tunable parameter. Then Neuberger's operator Z?n with mass m is given 
by 



m jj(Q) | m 



where 



L>W = (1 + s) U-A(A^A)- 1 



1/2 



(10) 



(11) 



What is important here in the definition of Neuberger's operator is the appearance of 
the square root of the operator A^A. This means that £>n connects all points of the lattice 
with each other. Note, however, that despite this the operator is a local operator in the 
field theoretical sense 15 . In practice, the operator A^A is represented by a matrix that 
is, unfortunately, very large. Having a physical volume of size L 4 , the number of sites is 
N 4 = (L/a) 4 . As the internal number of degrees of freedom per lattice point is 12, we 
end up with A being a (12iV 4 ) ® (127V 4 ) complex matrix with N typically in the range 
10 < N < 30 for present days simulations. Hence we have to construct the square root of a 
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Figure 1. The low-lying end of the eigenvalue spectrum of the matrix A. 

very large matrix. What is worse, to compute relevant physical observables, we don't need 
the operator itself but its inverse. Such an inverse is constructed by iterative methods 
like the conjugate gradient algorithm and its relatives 14 . 

The square root can be constructed by a polynomial expansion, normally based on a 
Chebyshev approximation, or by rational approximations. Both methods give comparable 
performances in practice. The convergence of the approximation to the square root is 
determined by the condition number of the (positive definite) matrix A^A. When the matrix 
A^A is normalized such that the largest eigenvalue is one, the condition number is given 
by the inverse of the lowest eigenvalue. We show in figure 1 the low-lying eigenvalues 
of A^A for different values of (3 where the parameter (3 is inversely proportional to the 
coupling strength of the theory. 

It can clearly be observed that very small eigenvalues can occur resulting in large con- 
dition numbers. In such situations the convergence of the approximations chosen can be 
rather slow and special tricks have to be implemented to accelerate the convergence. The 
most fruitful improvement is to treat a part of the low-lying end of the spectrum exactly 
by projecting this part out of the matrix A^A 16 > 17 . Further improvements can be imple- 
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mented, examples of which are discussed in ref. 17 . 

Despite all these technical improvements it is found that a typical value for the degree 
of a polynomial is 0(100) and a typical value for the number of iterations to compute the 
inverse of is again 0(100). Since in each iteration to compute D^ 1 the Chebyshev 
polynomial has to be evaluated, this means that for a value of a physical observable on 
a single configuration ten thousand applications of a huge matrix on a vector has to be 
performed. To compute the expectation value of some observable, this observables has to 
be averaged over many gluonic configurations. Clearly, this results in a very demanding 
computational effort and gives rise to a numerical challenge well suited for NIC. 

The only thing that helps a lot in this problem is that the matrix A is sparse. Only the 
diagonal and a few sub-diagonals are actually filled, a fact that finally makes the problem 
manageable - although, still a very large amount of computer resources are needed to 
tackle it. 

4 The scalar condensate 

The existence of an exact lattice chiral symmetry allows the use of finite size effects to test 
for SSB in QCD as in the case of the $ 4 -theory. Using a chiral invariant formulation of 
lattice QCD, it is possible to reach the region of very small quark masses where it is to be 
expected that chiral symmetry starts to get restored. 

The "magnetization" in the case of quantum chromodynamics is the condensate of a 
quark-antiquark state E = (V>V>). The role of the external magnetic field is played by the 
quark mass m q . We have developed a fully parallelized code with all technical improve- 
ments implemented. This allowed us to compute the scalar condensate as a function of 
the quark mass at several volumes. A (standard) caveat here is the fact that all compu- 
tations are done in the so-called quenched approximation where all internal quark loops 
are neglected. In QCD there is a peculiarity: the field configurations can have topolog- 
ical properties, characterized by the so-called topological charge which can be measured 
-unambiguously- through the number of zero modes of the operator D^. In fact, the 
formulae from quenched chiral perturbation theory are parametrized by the topological 
charge and it is hence very important to be able to identify the topological charge of the 
gauge field configurations. Without the special properties of lattice Dirac operators that 
satisfy the Ginsparg-Wilson relation such an identification would be very difficult. 

Let us give the complete theoretical formula from quenched chiral perturbation theory 
in lowest order: 

IU±i(m„ V) = S z [I v (z)K v {z) + I v+1 (z)K v -!{z)] + C ■ m q /a 2 . (12) 

The only important thing to notice here is that this relatively involved combination of 
Bessel functions do, as in the case of the $ 4 -theory, only depend on one scaling variable 

z = T,m q V (13) 

that contains the quantity of interest, namely the infinite volume, chiral limit scalar con- 
densate E. The additional term C ■ m q /a 2 is a power divergence that comes from the 
renormalization properties of the theory. We will not discuss this field theoretical aspect 
here but just notice that this term has to be included in the fit. 
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Figure 2. The scalar condensate computed on lattices of various size as a function of the quark mass. The solid 
lines are 2-parameter (S and the constant C of eq. (12)) fits according to the prediction of chiral perturbation 
theory. 



In figure 2, we show the result of our numerical computation of the scalar condensate 8 
in a fixed topological charge sector |^| = 1 as a function of the quark mass at several 
volumes. The solid line is a fit to the prediction of chiral perturbation theory, eq.(12). We 
find that the simulation data are described by this prediction very well. This means that 
we find evidence for the basic assumption on which the theoretical prediction relies: the 
appearance of spontaneous chiral symmetry breaking in (quenched) QCD. 

We want to remark that this work that has been performed at NIC was the first of this 
kind. The project consumed 1400 CPU hours on a typical distribution of the lattice on 128 
nodes. After this work, a number of other groups repeated such an analysis 21 > 22 ' 23 and 
it was reassuring to observe that very consistent results were found. In a subsequel work 
19,20 we developed a i so a quite general method for renormalizing the value of the bare 
scalar condensate as extracted from the finite size scaling analysis performed here. 
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5 Conclusion 



In this contribution we have demonstrated that by a combination of theoretical ideas, im- 
proved numerical methods and the use of powerful supercomputer platforms it is possible 
to test basic properties of field theories. Of particular interest was the question of whether 
the phenomenon of spontaneous symmetry breaking does occur in certain field theories 
important in elementary particle physics. The phenomenon of SSB leads to far reaching 
consequences in theories like QCD or the scalar sector of the electroweak interactions. In 
the work performed here, we found strong evidence for the appearance of spontaneous 
chiral symmetry breaking in quenched lattice QCD. This conclusion is the result of the 
fact that two theoretical concepts, the lattice approach to quantum field theories and chiral 
symmetry, met finally - and that enemies became friends. 

We are very much indebted to Hartmut Wittig for numerous discussions, suggestions 
and finally participating in the project at the stage when the non-perturbative renormal- 
ization of the scalar condensate was computed. L.L. thanks the INT at the University 
of Washington for its hospitality and the DOE for partial support during the completion 
of this work. This work was supported in part by the EU TMR program under contract 
FMRX-CT98-0169 and the EU HP program under contract HPRN-CT-2000-00145. 
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